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ABSTRACT 


semi -theoretical  model  of  the  upper  atmosphere  hao  beon  constructed 


to  allow  the  calculation  of  the  primary  atmospheric  properties  from 
90  to  10^000  km  altitude  as  a  function  of  local  time  and  solar  activity. 
The  calculations  are  based  upon  an  empirical  density  profile  from  a 
density  model  as  the  initial  input.  Assuming  diffusive  equilibrium 
above  110  km  and  isothermal  conditions  above  UOO  km,  the  remaining 
properties  (pressure,  temperature,  mean  molecular  mass,  and  the  individual 


particle  concentrations)  are  derived  by  using  the  hydrostatic  equation 


and  equation  of  state,  together  with  boundary  values  from  measurements. 

To  illustrate  the  model,  two  tables  are  given  that  present  the 
average  atmospheric  proper  Me  s  vs  altitude  for  sunspot  maximum  and  minimum 
conditions,  respectively,  for  the  months  immediately  following  the  Argus 

(September  1958)  and  Starfish  (July  and  August  1962)  high-altitude  nuclear 
detonations.  O  <■ 
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A  SEMI -THEORETICAL  MOIEL  FOR  ATMOSPHERIC  PROPERTIES 

FROM  90  to  10,000  KM 

I.  Introduction 

The  main  geophysical  problem  connected  with  atmospheric  structure 
in  the  hete  rosphere  (Table  1)  is  to  calculate  accurate  valu^-s  for  all  the 
primary  properties  as  a  function  of  time,  location,  and  date  by  taking 
into  account  the  relevant  processes.  The  primary  atmospheric  properties 
are  temperature,  density,  pressure,  and  mean  molecular  mass  (or  composition). 
To  account  for  their  variations,  the  following  factors  must  be  considered: 

(a)  time  (hour,  day,  sun-rotation  period,  season,  year,  sunspot  cycle); 

(b)  location  (altitude,  latitude,  longitude);  (c)  solar  activity 
(ultraviolet  radiation.  X-rays,  solar  plaama  and  associated  magnetic 
storms);  and  (d)  processes  (conduction,  diffusion,  mass  transport, 
photoionization,  dissociation,  recombination,  particle  escape  into  space). 

The  problem  of  describing  atmospheric  behavior  is  difficult  because  many 
of  the  above  elements  are  interrelated. 

Theoretical  attempts  have  been  made  to  form  a  model  of  the  high 
atmosphere,  that  is,  explain  the  observed  properties  in  terms  of  the 
processes  of  photoionization  and  recombination.  However,  these  fundamental 
attempts  have  been  vitiated  by  ignorance  of  the  solar  ionizing  spectrum 
and  uncertainties  in  the  absorption  and  ionizing  cross  sections  and  rate 
coefficients  for  the  varicis  possible  ionic  reactions.  For  example,  the 
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rate  coefficients  for  many  upper  atmosphere  reactions  are  not  known  by 
direct  measurement  but  are  inferred  from  certain  laboratory'  experiments 
conducted  at  various  equilibrium  conditions.  Consequently,  the  uncertainty 
in  the  rate  coefficients  extends  over  several  orders  of  magnitude.  For 
these  reasons,  present  data  and  theory  regarding  reactions  between  incoming 
solar  radiation  and  atoms  and  molecules  in  the  upper  atmosphere  are  not 
sufficient  to  provide  a  basis  for  accurate  calculations  of  atmospheric 
properties. 

The  magnitudes  of  various  atmospheric  properties  in  the  heterosphere 
can  be  derived  from  measurements  made  from  satellites,  rockets,  meteor 
observativons,  sky  emissions,  and  the  propagation  of  sound  and  radio  waves. 
Although  the  frequency  and  accuracy  of  the  various  measurements  are 
increasing,  the  data  still  are  few  and  contain  much  uncertainty.  The 
density  (drag)  data  resulting  from  tracking  satellites  are  generally 
considered  the  most  precise  and  are  by  far  the  most  numerous.  However, 
the  sparseness  of  the  data  in  time  and  space  requires  the  aid  of  a  model 
to  represent  the  variations  of  the  atmospheric  conditions  in  a  continuous 
fashion.  Therefore,  greater  understanding  of  the  upper  atmosphere  has 
been  attempted  through  an  accurate  description  of  all  phases  of  its 
complex  processes  by  means  of  models. 

A.  A  New  Atmospheric  Property  Model.  The  approach  followed  in 
the  construction  of  most  models  is  to  assume  the  altitude  dependence  of 
some  properties  in  order  to  calculate  those  remaining.  Most  of  the  models 
deal  principally  with  density,  pressure  and/or  temperature.  The  altitude 
variation  of  the  mean  molecular  mass  is  often  introduced  somewhat  arbitrarily, 
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and  therefore,  a  physically  consistent  vertical  distribution  of  the 

composition  cannot  be  obUined.  None  of  these  models  has  been  too  successful 
in  gaining  wide  acceptance. 

To  avoid  some  of  the  assumptions  of  other  methods  and  to  take  into 
account  some  factors  previously  neglected,  a  new  method  has  been  devised 
for  computing  atmospheric  properties.  In  this  work,  no  major  assumptions 
are  made  regarding  the  property  profiles.  Instead,  the.  primary  properties 
are  calculated  by  sUrting  with  an  empirical  density  profile  derived  from 
a  density  model  as  the  primary  input,  assuning  diffusive  equilibrium 
above  110  km  and  isotheimal  conditions  with  altitude  above  1,00  km.  The 
density  profile  is  represented  by  a  previous  empirical  model  (Anderson, 

1961»)  as  a  function  of  local  time  and  solar  activity  from  200  to  800  km. 

This  density  model  uses  corrected  values  of  the  10.7-cm  solar  flux 
(called  S  ')  as  an  index  of  the  sun's  extreme  ultraviolet  radiation  in 
order  to  account  for  the  marked  variations  in  density  that  occur  from 
day  to  day.  The  remaining  properties  are  derived  by  using  the  hydrosUtic 
equation  and  equation  of  state,  to  relate  density,  pressure,  temperature, 

and  mean  molecular  mass  to  altitude,  together  with  boundary  values  based 
on  measurements. 

In  summary,  to  develop  a  comprehensive  model  applicable  from  90  to 
10,000  km  the  theoretical  and  empirical  approaches  have  been  combined  in 
such  a  way  as  to  take  advantage  of  the  strong  points  of  each  by  maximizing 
the  number  of  constraints  on  the  atmospheric  properties  and  minimizing 
the  number  of  assumptions.  Since  the  calculations  are  quite  tedious, 
this  type  of  approach  is  feasible  only  with  the  use  of  modem  computer 
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techniques.  The  model  presented  in  this  report  has  been  programmed  for 
an  IBM- 7090  computer. 


IIj__Com£Osi^ion^_of__the_Jj£ger_^Atmosghere 

The  main  regions  of  the  earth's  atmosphere  are  presented  in  Table  1. 


TABU)  1 .  The  Main  Regions  of  the  Earth's  Atmosphere 


Atmospheric 

Region 

Dominant 

Process 

Altitude 
Range  (km) 

Remarks 

Homosphere 

Nixing 

Troposphere 

Stratosphere 

0-12 

12  -  50 

Mean  molecular 
weight  is 
constant 

Mesosphere 

50  -  90 

Heterosphere 

Diffusion 

Thermosphere 

90  -  550 

Frequent  particle 
collisions 

Exosphere 

550  -  60,000 

Collisions  rare; 
temperature  constant 
to  about  85CX)  km 

Since  the  dominant  process  which  controls  atmospheric  composition  up 
to  at  least  90  km  is  mixing,  the  composition  of  the  air  remains  the  same 
up  to  this  altitude  .^^He nee,  th '  mean  molecular  weight  (mass)  remains 
constant  up  to  90  km.  (Above  this  height,  the  molecular  weight  decreases 
as  the  composition  changes  with  altitude,  because  of  molecular  dissociation 
and  diffusion.)  Hence,  the  lower  altitude  boundary  for  the  model  was 
selected  at  90  km.  Below  this  altitude  the  U.S,  Standard  Atmosphere  (1962) 
is  perhaps  the  best  available  reference  for  atmospheric  properties.  A  mean 
molecular  weight  of  28.96  is  taken  from  this  reference  at  90  km.  This  value 
is  calculated  from  the  composition  of  dry  air  as  given  in  Table  2. 
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TABLE  2,  Normal  Composition  of  Clean,  Dry  Atmospheric  Air 
Near  Sea  Level  From  U.  S.  Standard  Atmosphere,  1962 


Constituent  Gas 

Content 

r  Cent  by  Volume) 

Molecular  Weight 

Molecular  Nitrogen 

78.0814 

28.0131* 

Molecular  Oxygen 

20,91*76 

31.9988 

Argon 

0.931* 

39.91*8 

Carbon  Dioxide 

0,0311* 

1*1*.  00995 

Neon 

0.001818 

20,183 

Helium 

0.000521* 

1*.026 

Krypton 

0.000111* 

83.80 

Xenon 

0,0000087 

131.30 

Molecular  Hydrogen 

0.00005 

2.01591* 

Methane 

0.0002 

16.01*303 

Nitrous  Oxide 

/ 

0.00005 

1*1*.0128 

Seasonal  and  latitudinal  variations  at  90  km  can  be  taken  into 
account,  if  desired,  by  using  values  from  the  supplemental  atmospheres 
derived  by  Court  et  al  (1962),  Table  3  gives  temperature  and  density  values 
for  the  U,  S,  Standard  Atmosphere,  1962, and  the  supplemental  atmospheres  at 
15°N,  30°N,  and  li5°N.  The  Tropical  atmosphere  represents  mean  annual  conditions 
at  15°N,  a  region  with  little  seasonal  variability;  the  Sub- tropical  and 
Mid-latitude  atmospheres  depict  mean  winter  and  summer  conditions  at  30°N 
and  li5°N,  respectively.  Since  the  U.  S.  Standard  Atmosphere,  1962,  is  a 
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middle«latitude  (approximate  US®)  year-round  mean  over  the  range  of  solar 
activity  between  sunspot  minima  and  maxima,  the  summer  and  winter  temperature 
and  density  values  given  for  the  Mid-latitude  (US^N  atmosphere  in  Table  3 

can  be  used  to  estimate  the  seasonal  variability  at  90  km  in  the  northern 

% 

hemisphere.  However,  for  the  model  given  in  this  report,  no  seasonal  or 
latitudinal  variability  is  taken  into  account  at  90  km.  Therefore,  a 
temperature  of  l8l®K  and  a  density  of  3.17  x  10"^  cm"^  are  taken  to 
apply  at  90  km. 


TABLE  3.  Atmospheric  Temperatures  and  Density  at  90  km 


Standard 

Temperature  (K) 

1 

Density  (am/cm  ) 

Atmosphere 

Summer  Winter 

Average 

Summer 

Winter 

Average 

U.S.  standard 
1962 

- 

l6l 

- 

- 

3.17x10"^ 

Hld'latltude 

(li5®N) 

17li  20li 

- 

3.60x10*^ 

3.32X10’*' 

- 

Sub- tropical 
(30°N) 

179  191 

- 

3.37x10"^ 

3.t»5xlO’^ 

- 

Tropical 

- 

I61t 

- 

- 

3.52x10"^ 

The  composition  of  the  upper  atmosphere  can  be  explained,  at  least  in 
a  qualitative  sense,  by  noting  the  types  of  photochemical  reactions  that 
can  occur.  In  the  folloviing  paragraphs,  the  reactions  leading  to  the  neutral 
constituents  will  be  discussed  first,  followed  by  the  ion  reactions. 
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A,  Oxygen  Dissociation  and  Recombination.  At  about  100  km,  the  absorption 
of  solar  radiation  with  wavelengths  shorter  than  l650  A  down  to  about  1300  A 
leads  to  the  dissociation  of  oxygen  molecules,  thus 

Op  ♦  hu  -  0  4  0  .  (1) 

For  wavelengths  shorter  than  1026  A,  the  oxygen  molecule  can  be  ionized, 

Og  4  hu  -  Og*  4  6.  (2) 

This  ionization  is  normally  followed  by  a  dissociative  recombination, 

Og*  +  e  -  0  +  0  .  (3) 

Although  the  oxygen  atoms  can  recombine  into  molecules,  photochemical 
equilibrium,  where  the  rates  of  dissociation  and  recombination  are  equal 
at  each  altitude,  does  not  prevail.  Instead  vertical  transport  due  to  both 
diffusion  and  mixing  plays  an  Important  role  in  determining  the  atomic  and 
molecular  concentrations  at  various  levels  near  100  km.  Nicole t  (1959)  has 
pointed  out  that  the  occurrence  of  oxygen  dissociatiai  in  a  region  where 
diffusive  transport  is  important  makes  it  impossible  to  deduce  either  the 
degree  of  dissociation  or  the  relative  concentration  of  atomic  oxygen  and 
molecular  nitrogen  from  measurements  of  molecular  oxygen  and  atmospheric 
density.  More  o-vgen  dissociates  than  recombines  above  ICX)  km,  due  to  the 
rapid  fall  off  with  altitude  of  the  recombination  processes.  Below  100  km 
collisions  occur  frequently  enough  for  recombination  to  prevail,  hence  more 
oxygen  recombines  than  dissociates.  Consequently,  there  is  a  steady  flux 
of  molecular  oxygen  upward  and  atomic  oxygen  downward  through  the  100  km 
level,  due  to  the  effects  of  diffusion  and  mixing.  Extreme  ultrav:  jlet 
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monochromator  measurements  by  Hinteregger  (1962)  show  that  atomic  oxygen 

becomes  the  dominant  atmospheric  constituent  above  approximately  150  tan 
altitude. 

B.  Nitrogen  Dissociation  and  Recombination.  The  most  active  process 
leading  to  the  dissociation  of  molecular  nitrogen  is  ionization  followed  by 
dissociative  recombination, 

N2^  +  e  -•  N  N  .  (1^) 

Atomic  nitrogen  can  react  with  molecular  oxygen  to  form  nitrous  oxide  and 
atomic  oxygen, 

N  Og  NO  +  0  . 

The  nitrous  oxide  in  turn  reacts  with  atomic  nitrogen  to  form  molecular 
nitrogen  and  atomic  oxygen, 

N  +  NO  -»  Ng  4  0  .  (5) 

The  effectiveness  of  these  reactions,  together  with  the  slowness  with  which 
molecular  nitrogen  dissociates,  causes  atmospheric  nitrogen  to  remain 
predominantly  in  molecular  form.  The  monochromator  measurements  by  Hinteregge 
(1962)  indicate  that  atomic  nitrogen  is  a  minor  atmospheric  constituent 
below  225-lan  altitude. 

C.  Atomic  Hydrogen.  Photodissociation  of  water  vapor  and  methane  near 

80  km  constitute  the  principal  sources  of  atomic  hydrogen.  Thus,  for  water 
vapor. 


HgO  +  hu  -  H  +  OH 


(7) 
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(8) 


This  reaction  is  usually  followed  by, 

OH  +  0  -  H  +  Og  . 

Ovdng  to  the  anaU  mass  of  the  hydrogen  atom  compared  to  other 

atmospheric  constituents,  the  hydrogen  concentration  does  not  decrease  with 

« 

altitude  as  rapidly  as  do  the  other  atmospheric  constituents  in  the  altitude 
region  where  diffusion  proceeds  rapidly;  hence  atomic  hydrogen  becomes  an 
increasingly  imporUnt  atmospheric  constituent  with  increasing  altitude. 
Atomic  hydrogen,  a  minor  constituent  in  the  thermosphere,  becomes  the 
dominant  constituent  above  1500  to  2500  km  altitude,  and  remains 
dominant  up  to  20,000  km  where  the  hydrogen  ion  becomes  dominant. 

The  distribution  of  atomic  hydrogen  out  to  60,000  km  from  the  earth 
at  the  extremes  of  the  sunspot  cycle,  according  to  Johnson  (1961a)  is 
shown  in  Fig.  1.  Two  curves  are  shown,  one  for  the  maximum  of  the  simspot 
cycle  (average  temperature  of  1500  K)  and  one  for  the  minimum  (average 
temperature  of  1000  K).  The  source  of  atomic  hydrogen  near  80  km  can  be 
expected  to  remain  essentially  constant  through  the  sunspot  cycle,  but  the 
rate  of  escape,  depending  on  the  temperature  at  the  base  of  the  exosphere 
varies  with  the  sunspot  cycle.  The  escape  will  be  relatively  rapid  when 
the  temperature  is  high,  and  the  concentration  of  hydrogen  will  be 
correspondingly  low  in  the  exosphere  near  sunspot  maximum.  The  escape  is 
relatively  slow  when  the  temperature  is  low,  so  that  the  concentration  must 
be  comparatively  high  near  sunspot  minimum. 

D.  telium.  Nicolet  (1961)  showed  that  helium  atoms  are  an  important 
constituent  in  the  lower  exosphere.  He  explains  the  high  densities  derived 
from  the  rate  of  change  of  the  period  of  Echo  satellite  by  the  presence  of 
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helium.  Evidently,  atomic  oxygen,  nitrogen  or  hydrogen  cannot  explain 
the  slow  density  decrease  between  750  and  1500  km.  Although  atomic  oxygen 
is  the  most  important  constituent  in  the  upper  thermosphere,  atomic  helium 
dominates  over  atomic  oxygen  somewhere  above  800  to  1000  km.  Atomic  hydroge 
dominates  over  helium  somewhere  above  1500  to  2500  km. 

E.  Ionospheric  Composition.  The  primary  ions  formed  in  the  E  region 
(85  to  lliO  km)  of  the  ionosphere  are  02*^,  and  o"^.  The  N2*^  ions 
dissociatively  combine  very  rapidly,  and  they  may  also  react  with  oxygen, 

so  that  the  concentration  of  N2  is  small.  Ion-atom  exchange  reactions  of 
the  type 

yz  *  XY*^  +  Z  (9) 

proceed  very  rapidly,  so  that  the  reactions  of  0^  with  0^  and  N2  quickly 
remove  the  o'*’  ions  and  produce  either  NO^  or  O2*.  The  dominant  ions  in 
the  region  (11*0  to  200  km)  are  NO*  and  Og*  near  the  lower  altitude 

boundary,  with  a  gradual  transition  to  0*  as  the  principal  ion  at  the  upper 
boundary.  0*  becomes  the  dominant  ion  because  of  the  rapid  decrease  with 
altitude  of  the  neutral  molecular  constituents  that  otherwise  would  tend 
to  eliminate  0*  by  reaction  (9).  In  the  F2  region  (200  to  8OO  km),  the 
ions  present  are  0  and  N  ,  with  0*  greatly  predominant.  The  helium  ion 
starts  to  become  dominant  between  8OO  and  ll*00  km  (Bourdeau  and  Bauer, 

1963),  depending  on  atmospheric  temperature.  Likewise,  protons  become  the 
dominant  ion  between  about  11*00  to  1*000  km.  Thus,  there  are  two  transition 
regions  (from  oxygen  to  helium  ions  and  from  helium  to  hydrogen  ions)  in 
the  upper  ionosphere  rather  than  a  single  transition  from  oxygen  to  hydrogen 
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as  previously  believed.  The  outer  portion  of  the  earth's  atmosphere 
(extending  from  about  20,000  km  to  at  least  U0,000  km  above  the  earth's 
surface)  consists  mainly  of  protons.  The  protons  predominate  in  the  earth's 
outer  atmosphere  over  the  hydrogen  atoms  because  the  confining  effect  of 
the  earth's  magnetic  field  makes  their  rate  of  escape  to  space  much  less  than 
for  the  neutral  hydrogen  atoms. 

III.  The  Physics  of  the  Upper  Atmosphere 

Assuming  that  the  earth's  atmosphere  is  a  continuous  medium  consisting 
of  a  gas  in  static  equilibrium,  the  equation  connecting  the  pressure  P  and 
density  p  at  any  altitude  h  is 

dP  .  -g  f  dh  ,  (10) 

where  g  is  the  gravitational  acceleration.  The  variation  in  g  with  altitude 
can  be  taken  into  account  by  using  the  relationship 

g  -  go  RgARo  +  h)^,  (11) 

where  g  is  the  acceleration  of  gravity  at  the  earth's  surface  and  R  is 

o  o 

the  earth's  radius.  Eq.  10  is  the  hydrostatic  equation  in  differential 
form.  Further,  let  m  be  the  mean  iiolecular  mass  of  the  gas  and  n  denote 
the  number  density  or  concentration;  the  density  p  is  then 

p  =  n  m  .  (12) 

Since  the  terrestrial  atmosphere  is  almost  a  perfect  gas,  the 
pressure  is  given  by  the  equation  of  state 

P  =  n  k  T  ,  (13) 
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where  T  is  the  temperature  and  k  the  Boltzmann  constant  (1.30  x  lO"^^  erg 
<l8g'^). 

Multiplying  both  sides  of  (12)  by  g  and  dividing  the  result  into  (13), 

gives 

P/p  g  *  kT/mg  .  H,  (11^) 

where  H  -  kT/mg  is  known  as  the  pressure  scale  height,  a  parameter 
convenient  for  atmospheric  calculations.  It  should  not  be  confused  with 
the  density  scale  height,  -P/(dp/dh),  a  parameter  used  to  relate  density 
to  satellite  drag  (Anderson,  1961*). 

Above  liOO  km  ions  are  distributed  independent  of  the  neutral  particles 
as  a  result  of  diffusive  equilibrium.  The  scale  height  of  a  plasma 
(an  electrically-neutral  combination  of  ions  and  electrons)  containing 
a  single  ionic  species  of  mass  m^  and  temperature  T^^  is,  neglecting  the 
mass  of  the  electron  with  respect  to  the  ions,  from  (lli). 

Hi  .  k  (Ti  ♦  Tg)/  mig  ,  (15) 

where  T^  is  the  temperature  of  the  electrons. 

Even  though  immediately  after  formation  the  ions  may  have  more  than 
thermal  energy,  they  speedily  lose  the  excess  by  colliding  with  neutral 
particles,  so  that  their  temperature  becomes  the  same  as  that  of  the  neutral 
particles.  This  may  not  be  true  for  the  electrons  because  of  the  large 
mass  difference.  However,  if  the  ions  and  electrons  are  in  thermal 
equilibrium,  T^  «  T^  «  T  in  (15)  and  the  ions  and  neutrals  have  the  same 
mass,  the  scale  height  for  the  ions  and  electrons  is  then 
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•  2  kT/ng  . 


(16) 


Thus,  the  ion  scale  height  is  twice  that  for  neutral  particles  of  the  same 
mass.  The  electric  field  established  between  the  electrons  and  ions  just 
cancels  half  the  gravitational  force  acUng  on  the  ions.  When  there  are 
several  kinds  of  ions,  the  average  electric  field  is  determined  by  the  mean 
ionic  mass  m.  For  an  ion  of  mass  m,  the  effect  of  the  electrostatic  field 
required  to  maintain  charge  neutrality  is  equivalent  to  chan;;ing  gravity  by 
the  factor  (1  -  in/Zm).  All  ionic  constituents  will  be  subjected  to  the  same 
average  electric  field.  Mange  (I960)  gives  a  method  for  calculating  the 
altitude  distribution  of  ionic  constituents  in  electrostatic  equilibrium. 

Eq.  10  to  16  are  applicable  in  the  atmosphere  where:  (a)  the  escape 
into  space  of  a  constituent  is  not  important,  (b)  the  earth's  magnetic  field 
does  not  change  the  vertical  distribution,  and  (c)  the  effect  of  the  earth's 

axial  rotation  is  negligible.  Each  of  the  above  factors  will  now  be 
consi(jered  in  turn. 

— Accuracy  of  the  Hydrostatic  Assumption.  The  hydrostatic  equation  (lO) 
IS  based  on  the  concept  of  local  thermodynamic  equilib.'ium  in  which  collisions 
between  particles  are  frequent  enough  so  that  there  is  a  complete  statistical 
exchange  of  particle  energy  and  momentum  in  the  volume  under  consideration. 

In  the  thermosphere  the  particles  collide  often  enough  to  insure  a  Maxwellian 
velocity  distribution  and  the  existence  of  a  meaningful  kinetic  temperature. 

In  the  exosphere,  above  about  500  km,  collisions  are  sufficiently  reduced 
so  that  the  existence  of  a  Maxwellian  distribution,  and  consequently,  the 
validity  of  applying  the  hydrostatic  relationship,  must  be  examined  anew. 
Liouville's  theorem  can  be  used  to  show  that  the  Maxwellian  velocity 
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distribution  prosent  bolow  the  b&se  of  the  exosphere  Applies  t'Qually  well 
in  the  exosphere,  provided  that  the  escape  of  particles  to  space  is  negligible 
(Herring  and  Kyle,  1961),  The  atmosphere  is  so  rarefied  in  the  exosphere 
that  there  is  very  little  solar  railiation  absorbed.  Also,  there  are  no 
energy  loss  mechanisms  to  disturb  the  Maxwellian  velocity  distribution 
(Sharp,  1962).  If  the  neutral  particles  have  a  Maxwellian  velocity 
distribution,  the  ion  velocities  must  also  have  the  same  distribution  and  a 
meaningful  kinetic  temperature.  Thus,  if  the  ion  and  neutral  particle 
temperatures  are  the  same  in  the  thermosphere,  they  must  also  be  the  same 
in  the  exosphere  wherever  the  neutral  particles  have  a  Maxwellian  velocity 
distribution.  At  the  temperatures  present,  the  distribution  is  disturbed 
only  by  the  escape  of  hydrogen  from  well  above  the  base  of  the  exosphere. 

Since  at  very  high  altitudes,  (a)  the  hydrogen  ion  density  is  comparable 
to  or  higher  than  the  neutral  hydrogen  density,  and  (b)  the  collision  cross 
section  for  the  ions  is  much  larger  than  that  for  the  neutral  particles,  it 
follows  that  the  maximum  altitude  in  the  atmosphere  where  the  hydrostatic 
equation  is  applicable  can  be  determined  by  calculating  the  altitude  in  the 
exosphere  where  the  density  of  the  neutral  hydrogen  starts  to  be  significantly 
affected  by  the  escape  of  hydrogen  atoms, 

Johnson  (1961b)  divides  the  Maxwellian  velocity  distribution  of 
exospheric  particles  into  four  groups  depending  on  their  trajectories  with 
reference  to  the  base  of  the  exosphere.  These  four  groups  are;  (1)  particles 
with  less  than  escape  velocity  that  have  come  up  to  the  level  concerned  on 
elliptic  paths  from  the  base  of  the  exosphere;  (2)  particles  with  less  than 
escape  velocity  whose  perigees  lie  entirely  above  the  base  of  the  exosphere 
(these  particles  are  in  trapped  orbits);  (3)  particles  with  more  than  escape 
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volocity  whose  perigees  lie  below  the  base  of  the  exosphere  (half  of  these 
particles,  designated  group  3a,  are  in  the  process  of  escaping,  having 
attained  escape  velocity  in  a  collision  near  the  base  of  the  exosphere; 

Uio  remaining  half,  designated  group  3b,  are  in  the  process  of  being 
captured  from  space);  (I4)  particles  with  more  than  escape  velocity  whose 
f)erigees  lie  above  the  base  of  the  exosphere  (these  are  particles  from 
.jpaci?  that  approach  the  earth  less  closely  than  the  base  of  Uie  exosphere). 
The  corresponding  r  )rmalized  groups  are  designated  Xj,  X2,  Ky  and 
r.'jpecti  vely.  On-  of  us  has  expressed  these  fractional  groups  in  terms  of 
error  functions  in  Appendix  A  of  this  report.  These  analytical  expressions 
for  the  above  distributions  (Eq.  1,  2,  3,  and  k  in  Appendix  A)  can  be  used 
to  calculate  exospheric  densities  out  to  about  5  earth  radii.  Fig.  2 
presents  curves  derived  from  the  results  of  calculations  made  for  the 
fraction  of  particles  present  in  a  Maxwellian  distribution  and  also  present 
in  the  atmosphere  as  a  function  of  altitude  h  above  the  earth's  surface, 
out  to  10,000  km  for  three  selected  temperatures,  875K,  1125k,  and  1625k. 

The  curves  start  at  h  =  hOO  km  since,  for  reasons  given  later,  the  atmospheri 
temp,  nature  above  1^00  km  over  ,  given  location  on  the  earth's  surface  does 
not  vary  wi  tii  altitude.  Consequently,  for  a  given  solar  activity,  the  curves 
could  be  labeled  also  according  to  local  times  instead  of  temperatures.  The 
tempcratur-'  of  875K  corresponds  approximately  to  a  nighttime  temperature 
near  sunspot  minimum  conditions,  while  the  temperature  of  1625K  is  a  day 
time  temperature  during  sunspot  maximum  conditions  (Anderson,  19614).  The 
1129k  temperature  represents  a  2i4-hr.  mean  for  average  solar  activity 
Condi tions. 


IS 
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Two  sets  of  curves  are  shown  in  Fig.  2  for  the  two  cases:  (a)  l-x^> 

-y^/2,  since  in  the  actual  distribution  groups  3b  and  i*  are  clearly  absent, 
as  the  hydrogen  concentrations  in  interplanetary  space  are  too  small  to 
contribute  significantly  to  the  hydrogen  captured  by  the  earth  and  (b) 

“X3/2  -X2>  wiih  particles  in  group  2  absent.  With  regard  to  group  2 
particles,  Brandt  and  Chamberlain  (I960)  argue  that  they  are  present  (case  a), 
wfiile  Opik  and  Singer  (I960)  argue  that  they  are  not  (case  b).  Since  the 
results  of  the  calculations  made  for  the  two  cases  are  markedly  different, 
both  cases  have  been  shown  in  Fig.  2  for  comparison.  If  it  is  assumed  that 
the  hydrostatic  equation  can  no  longer  be  used  to  make  accurate  calculations 
for  those  conditions  in  which  the  fraction  of  Max'iellian  particles  is 
below  0.90,  then  for  an  average  b^mperature  of  1125K  it  can  be  applied 
accurately  out  to  8^00  km  if  case  (a)  holds  or  only  out  to  2900  km  if  case 
(b)  holds.  However,  Johnson  (1961b)  has  considered  both  cases  and  suggests 
it  is  likely  that  the  trapped  orbits  are  populated  out  to  6  or  7  earth 
radii.  Since  some  collisions  occur  above  the  base  of  the  exosphere,  there 
will  be  some  injection  of  particles  into  trapped  orbits  in  the  lower  exosphere. 
Also,  the  presence  of  ionized  hydrogen  with  a  Maxwellian  velocity  distribution 
out  to  5  or  6  earth  radii  (Johnson,  1960b)  contributes  to  the  population  of 
trapped  orbits  for  neutral  particles;  charge  exchange  can  occur  between  the 
neutral  hydrogen  (principally  group  1  particles)  and  the  ionized  hydrogen, 
tending  to  fill  the  trapped  orbits.  Therefore,  following  Johnson,  the 
constraints  imposed  by  case  (a),  rather  than  case  (b),  were  adopted  for 
tho  model.  Fig.  2  indicates  that  the  error  involved  for  either  case  will 
increase  with  altitude  and  temperature.  The  greatest  error  in  using  the 
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hyfiroutatic  a.'j^JunipMon  will  occur,  as  illustrated  in  Fig. 2,  during  the  day 
near  sunspot  maximuin  conditions.  For  case  (a)  the  error  will  start  to 
hf'como  significant  above  about  ^000  km  above  the  earth's  surface, 

B.  Kffo‘ct  fjf  tiie  Earth's  Magnetic  Field  on  the  Vertical  Distribution 
ol‘  F'artlcles.  Th';  magnetic  field  affects  the  distribution  of  ions,  but  not 
distribution  of  neutral  particles.  .Since  the  mean  free  jjath  for  tho 
f)rotons  ar  i  (,‘loctrons  ar*.-  long  compared  with  the  cyclotron  radii,  the 
aV(>r’ago  motion  of  the  protons  and  electrons  is  along  the  magru?ti.c  field 
liri'  and  tli.  rn'.an  motion  of  the  particles  is  not  isotropic.  However, 
.Johnson  (lV60b)  points  out  that  the  detailed  velocity  distribution  remains 
i;:otro[jic,  oven  if  the  mean  motion  of  the  particles  is  not,  if  collisions 
■among  the.'  prtjtons  and  eJectrons  remain  rather  frequent.  Out  to  10,000  km, 
at  lc>ast,  collisions  do  remain  frequent  enough,  on  the  order  of  one  per 
.ocond,  ev('ri  i’or  particie  conccnitrations  and  night  time  bjmpera turns  prf'serit 
n-ar  solar  minimum.  The  distribution  of  ions  must  follow  Uu;  hydrostatic 
rvdati  onshif)  so  long  as  the  detailc.Mj  velocity  distribution  remains  isotropic 
(;'pi  t//,''r,  IVS'2). 


O.hffcjcts  (jf  the  Farth'.s  Axial  Rotation  on  Ion  Distribution.  The  magn<-Lic 
fi'  ld  aff(,'cls  the  di.itribution  of  ions  by  linking  the  if)n-olec tron  pla;.7n.i 
in  thf  '.‘XcjSfjlu.'ro  to  the'  iowt'r  ionosphere,  thus  requiring  that  the  pl;.sma 
r(jtat(;  with  tln'  earth.  The  resulting  centrifugal  force  on  the  ions  modifies- 
the  gravitational  attrat  tion  that  enters  into  the  hydrostatic  equation. 

.-inc;  all  th"  ions  and  electrons,  wherever  they  may  be  in  tlu'  magn(?tic 
field,  must  rotate  with  the  earth,  the  generalized  hydrostatic  equation 
must  be  corrected  for  the  earth's  rotation.  The  modifio'd  hydrostatic  equ.-ition 
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for  tho  distrib'jtion  of  ions  or  electrons  in  a  region  where  diffusive 
equilibrium  prevails  is  (Johnson,  1960b) 


N(h)  =  exp  m/(R^+h)  2kTL  exp  tm  fj  ^  cos^fl/likT],  (17) 

where  N(h)  is  th-  number  den.iity  at  altitude  h,  is  a  normalization 
constant,  g  is  the  acceleration  of  gravity  at  the  earth's  surface,  R 

o 

IS  the  radius  of  the  earth,  m  is  the  ion  mass,  n  is  the  rotational  velocity 
of  the  earth  about  its  axis,  e  is  the  geocentric  angle  between  the  equatorial 
plane  and  the  point  considered  (i.e.  latitude),  and  Ic  the  Boltzmann  constant. 
Th=?  first  exponential  is  the  normal  relationship,  and  the  second  exponential 
represents  the  modification  due  to  the  centrifugal  force  (  At  h  =  10,000  km. 


the  maximum  value  of  the  second  exponential  occurs  over  the  earth's  equator 
and  is  inversely  p»’oportional  to  the  temperature.  For  T  =  800K,  it  is  1.0^^# 
D.  Diffusive  Equilibrium.  Inasmuch  as  the  atmosphere  consists  of  a 
mixture  of  gases  subject  to  a  force  field,  the  equilibrium  distribution  of 
its  constituents  can  be  expected  to  show  some  degree  of  diffusive  separation. 
Although  turbulent  mixing  below  100  km  does  not  allow  the  development  of 
diffusive  equilibrium  distributions,  above  this  altitude  there  is  experimental 
evidence  that  diffusive  separation  commences.  In  the  geopotential  field, 


diffusive  equilibrium  gives  a  concentration  distribution  for  each  neutral 
species  that  varies  exponentially  with  the  geopotential,  with  a  more  rapid 
decrease  in  concentration  with  increasing  geopotential  for  the  heavier 
constituents  in  the  atmosphere  than  for  the  lighter  (Chapman,  I960).  The 
result  is  a  static  distribution  of  gas  constituents  under  the  action  of 
the  external  force  field,  expressed  by  the  hydrostatic  relationship  applying 
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for  each  atmospheric  constituent  independently  of  the  others.  This 
generalization  holds  for  the  isothermal  conditions  that  always  exist  above 
1|00  km.  Below  this  altitude,  the  temperature  gradient  present  introduces 
the  effect  of  thermal  diffusion,  which  slightly  modifies  the  conclusion. 
From  (10),  (12),  and  (13), 

-Pi 

dP^  «=  -  g  dh  -  -n^  m^  g  dh  =  -  m^  g  dh  ,  (18) 

kT 

where  the  subscript  1  refers  to  a  particular  neutral  constituent  having 
density  concentration  n^,  molecular  mass  m^,  and  partial  pressure 
Likewise,  for  another  constituent  whose  concentration  is  n^,  molecular  mass 
m2,  and  partial  pressure  P2. 

^  =  -  gdh  (19) 

m2P2  kt 


If  both  constituents  are  at  the  same  altitude  h,  then 


dP, 


m2P2 


(20) 


Integrating  both  sider.  of  (20),  the  left  side  from  pressure  P2(h^)  to 
pressure  P2(h2),  corresponding  to  altitude  h^  and  h2,  respectively,  and 
the  right  side  from  Pi(hj^)  to  P^(h2),  gives 


^2^^2^dP« 
o  c 


P2(hi) 


n,  dp 

C  f\  1 


m. 


or 


Pl(hi) 
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Thus,  if  the  pressure  altitude  dependence  of  one  neutral  constituent  is  known 
the  others  can  be  detennined  by  (21).  If  j  is  the  total  number  of  consti¬ 
tuents,  it  follows  from  Dalton's  law  that 

>  (22) 

1 
and 

J 

y  -  P(h2)  ,  (23) 

lAiere  P(h^)  and  P(h2)  are  the  total  pressures  at  altitudes  h^  and  h2, 
respectively.  If  the  partial  pressures  are  known  for  each  of  the  J 
constituents  at  one  pressure  altitude  P(hj^),  they  can  be  determined  for  any 
other  pressure  altitude  P(h2)  from  (21)  and  (23)  provided  that  diffusive 
equilibrium  exists  from  h^  to  h2. 

With  regard  to  ions,  Johnson  (1960b)  points  out  that  above  the  Fg 
maximum  in  the  ionosphere  (300-350  km),  diffusion  becomes  so  rapid  that  it, 
rather  than  recombination  becomes  the  dominant  factor  controlling  the  ion 
dis'c  )Ution.  Diffusion  of  ions  through  the  neutral  gas  tends  to  distribute 
the  ions  so  that  they  satisfy  a  hydrostatic  law  (see  paragraph  following 
Eq.  16)  Independently  of  the  neutral  particles. 
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During  1956-1958  the  Naval  Research  Laboratory  launched  four  rockets 
instrumented  with  radio-frequency  mass  spectrometers  to  study  neutral  gas 
compositijn  above  100  km  at  Churchill,  Canada  (59°N).  Results  from  the 
flights,  two  at  night  and  two  during  the  day,  showed  that  diffusive  separation 
started  in  the  approximate  range  from  100  to  120  km  (Meadows  and  Townsend, 
I960).  USSR  measurements  made  during  the  night  in  September  I960  indicated 
that  the  height  of  the  beginning  of  diffusion  was  about  110  km  (Pokhunkov, 
1962).  In  view  of  the  foregoing,  the  altitude  where  diffusive  equilibrium 
begins  is  tentatively  taken  at  110  km  for  the  model  in  this  report.  This 
altitude  may  be  changed  later  when  more  knowledge  becomes  available 
concerning  diurnal,  latitude,  seasonal,  and  solar  effects  on  the  exact 
altitude  where  diffusive  separation  starts. 

E.  Isothermal  Conditions.  The  concept  that  the  atmosphere  extmdlng 
above  1^00  km  over  a  given  location  on  the  earth's  surface  is  isothermal 
in  the  sense  that  the  temperature  does  not  vary  with  altitude  is  now  well 
established.  In  the  tenuous  gas  of  the  upper  thermosphere  the  thennal 
conductivity  is  independent  of  the  pressure  while  the  heat  capacity  varies 
linearly  with  density.  Consequently,  the  conductivity  is  very  large  compared 
to  the  heat  capacity  (Johnson,  1956).  Above  hOO  km,  the  absorption  of 
energy  is  negligible  and  the  relatively  high  heat  conductivity  eliminates 
temperature  differences;  hence  the  temperature  is  nearly  constant  with 
altitude  for  many  thousands  of  kilometers.  Kicolet  (I960)  has  concluded 
that  the  time  of  conduction  is  so  short  that  there  is  a  strong  tendency  to 
a  vertical  isothermy  for  the  atmosphere  at  high  altitudes  above  a  given 
region  of  the  earth. 
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Spitzer  (19ii9)  first  pointed  out  that  the  outermost  atmosphere  or 
exosphere,  the  region  in  the  atmosphere  in  which  collisions  are  relatively 
rare,  should  be  isothermal,  that  is,  that  the  distribution  of  atoms  in  the 
exosphere  should  be  that  of  an  isothermal  gas.  The  base  of  the  exosphere  is 
defined  as  the  level  at  which  the  mean  free  path  is  equal  to  a  scale  height 
H  (Eq.  lU).  Its  altitude  usually  lies  between  hSO  to  650  km.  The  upper 
extent  of  the  exosphere  may  be  taken  to  coincide  with  the  outer  edge  of  the 
magnetosphere,  extending  out  to  at  least  U0,000  km  above  the  earth's  surface 
(Beard,  I960). 

The  principal  heat  source  in  the  upper  atmosphere  above  130  km  is  the 
extreme  ultraviolet  radiation  (EUV)  from  the  sun  having  wavelengths  from 
2CX)  to  900  A.  Magnetic  (corpuscular)  heating  is  not  important  under  normal 
conditions,  but  during  the  disturbed  conditions  accompanying  major  magnetic 
storms,  it  can  lead  to  transient  heating  increases.  Fig.  3  (Friedman,  1959) 
shows  the  altitude  at  which  a  fraction  1/e  remains  of  the  radiation  incident 
on  the  atmosphere.  Almost  all  of  the  energy  of  the  incident  EUV  radiation 
is  absorbed  by  the  main  atmospheric  constituents  Og#  0  and  Ng  in  the  altitude 
region  from  100  to  200  km.  The  radiation  absorbed  below  200  km  is  responsible 
for  heating  the  atmosphere  at  altitudes  extending  many  thousands  of  kilometers 
above  200  km  (Johnson,  1958) •  The  degree  of  heating  (represented  by  the 
temperature  of  the  isothermal  atmosphere)  is  proportional  to  the  energy 
absorbed  in  the  region  between  100  and  200  km.  The  amount  of  energy 
absorbed  depends  on  its  intensity  near  the  top  of  the  absorbing  region, 
which  is  a  function  of  (a)  solar  activity,  and  (b)  the  location  of  the 
absorbing  region  with  respect  to  the  sun-earth  line.  Hence,  the  portions  of 

22 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


the  exosphere  over  different  regions  of  the  earth  each  will  have  different 
isothermal  temperatures.  Anderson  (1961*)  takes  (a)  into  account  by  using 
the  parameter  S'  as  an  index  of  the  EUV  emitted  from  the  sunj  the  effect 
of  (b)  is  described  by  using  t,  the  local  time,  as  a  parameter.  These  two 
parameters  can  be  used  also  to  represent  the  variation  of  the  exospheric 
temperature  in  the  night  hemisphere,  inasmuch  as  the  nighttime  temperature 
in  a  given  region  of  the  exosphere  depends  on  the  temperature  existing 
before  sunset  in  the  corresponding  lower  absorbing  region  of  the  thermosphere 
(a  function  of  S')  and  on  the  rate  of  cooling  of  this  lower  region  with 
time  after  sunset  (a  function  of  t).  An  estimate  of  the  temperature  in  the 

upper  thennosphere  and  lower  exosphere  as  a  function  of  S'  and  t  is  shown 
in  Fig.  1*. 

The  altitude  of  the  thennopause  is  subject  to  a  considerable  diurnal 
variation  varying  with  the  amount  of  EUV  abosrbed  between  100  and  200  km;  it 
is  a  maximum  in  the  sunlit  atmosphere  and  a  minimum  in  the  dark  atmosphere. 
During  the  day  at  sunspot  maximum,  the  altitude  of  the  thermopause  is  about 
liOO  km;  during  the  night  at  sunspot  minimum,  its  altitude  is  about  200  km. 
Most  of  the  solar  energy  absorbed  below  200  km  is  conducted  downward  to 
below  120  km  where  it  is  dissipated  by  infrared  radiation.  Thus,  there  is 
a  steep  positive  temperature  gradient  between  100  and  200  km  ending  at  the 
thermopause  or  the  level  above  which  the  atmosphere  becomes  isothermal. 

Above  120  km,  the  atmosphere  is  a  very  poor  radiator  in  the  infrared,  but 
below  100  km,  the  atmosphere  is  a  much  better  radiator  becfuse  of  the  presence 
of  H2O,  OH,  COg  and  0^. 
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The  isothermal  temperature  used  in  the  model  is  the  kinetic  temperature 
that  can  be  detennined  only  for  a  gas  with  a  Maxwellian  velocity  distribution. 
As  pointed  out  under  the  discussion  of  the  accuracy  of  the  hydrostatic 
assumption)  the  Maxwellian  velocity  distribution  applies  in  the  exosphere 
provided  that  the  escape  of  particles  to  space  is  negligible.  For  hydrogen, 
the  escape  of  atoms  is  comparatively  rapid  so  that  the  velocity  distribution 
of  hydrogen  atoms  in  the  upper  exosphere  is  not  Maxwellian.  Therefore,  the 
hydrogen  atoms  in  the  upper  exosphere  have  a  non-Maxwellian  distribution 
that  becomes  more  pronounced  with  altitude.  Under  these  circumstances  the 
concept  of  kinetic  ten5>erature  is  not  entirely  applicable,  although  an 
effective  temperature  can  be  defined  by  considering  the  average  energy  of 
the  hydrogen  atoms.  This  effective  temperature  is  not  constant  as  a  function 
of  distance  from  the  earth.  The  same  altitude  above  which  the  use  of  the 
hydrostatic  equation  becomes  questionable  in  the  model  can  also  be  taken  to 
be  where  the  deviation  between  the  kinetic  and  effective  temperature 
becomes  significant.  This  altitude,  above  »4iich  the  fraction  of  Maxwellian 
particles  was  below  0.90,  was  at  8500  km  for  an  average  temperature  of  1125  K, 
assuming  particles  are  present  in  trapped  orbit.  Sharp  (1962,  Fig.  5) 
indicates  the  effective  temperature  decreases  with  distance  out  to  10  earth 
radii  for  an  assumed  kinetic  temperature  of  1000  K  at  the  base  of  the  exosphere. 

JV,  The  Impirloai  flsnaity  Model 

In  Mils  imuMon,  the  empirical  density  model  used  to  derive  the  primary 
ftir  the  simitepherlo  property  model  will  be  described  in  brief.  It 
h  «eiiteil  In  ilntail  elsewhere  (Anderson,  1962,  1961*).  The  density  p, 

I  '  \  In  repreiNinied  by 

21* 


IMI 


LOC 


CK  COMPANY 


(21*) 


f  -  K"  S  , 

where  K»'  and  N  are  constants  for  a  given  altitude  h  and  local  time  t,  and 

S'  is  a  daily  index  of  the  sun's  extreme  ultraviolet  radiation  expressed 

“22  “2  “1 

in  units  of  10  watt  m  cps  ,  K"  and  N  are  presented  in  tables 
(Anderson,  196U)  for  selected  altitudes  from  200  to  800  km  and  every  hour. 
Hence,  using  (21^)  together  with  the  tables,  the  density  and  its  variations 
can  be  readily  calculated  as  a  function  of  S',  derived  from  the  relative 
variations  of  the  sun's  10.7-cm  flux  (S)  and  satellite  drag  data.  Table  h 
gives  average  values  for  S'  and  S  by  months,  1958  to  1962.  Eq.  2h  has  been 
programmed  for  an  IBM-7090  computer;  the  K”  and  N  values  corresponding  to 
a  given  h  (nearest  kilometer)  and  t  (nearest  h  minutes)  are  found  from  the 
tables  by  logrithmic  interpolation  for  K”  and  linear  interpolaUon  for  N. 
Since  the  results  of  previous  work  indicate  that  density  variations  with 

latitude  are  anall,  at  least  between  60  N  and  60  S,  they  have  been  neglected 
in  this  model. 
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TABLE  li. 


Average  S' 

and  S  (10 

-22 

-2 

watt  m 

cps  ^ 

)  by  months. 

1956  to  1962. 

Month 

1956 

1959 

I960 

1961 

1962 

S' 

s 

S' 

S 

S' 

S 

S' 

s 

S' 

s 

January 

- 

- 

2LS 

271 

199 

200 

no 

120 

55 

93 

February 

251» 

210 

252 

206 

I8l 

169 

9a 

105 

62 

101 

March 

279 

250 

256 

226 

17a 

146 

no 

104 

73 

100 

April 

322 

246 

2ia 

210 

169 

167 

92 

105 

81 

96 

May 

272 

219 

209 

213 

1^7 

163 

81 

99 

62 

98 

June 

206 

220 

m 

217 

133 

162 

75 

no 

55 

91 

July 

212 

224 

160 

203 

130 

164 

78 

116 

a2 

81 

August 

221t 

237 

167 

234 

153 

17a 

78 

106 

a6 

77 

September 

272 

243 

207 

194 

169 

164 

82 

112 

65 

89 

October 

266 

226 

169 

164 

163 

141 

79 

96 

70 

87 

November 

2ii9 

207 

209 

163 

160 

ia7 

69 

69 

65 

8a 

December 

265 

234 

169 

179 

13a 

136 

65 

93 

56 

81 

All  Months 

258 

226 

212 

208 

161 

161 

8a 

104 

61 

90 

V.  The  Model  for  Atmospheric  Properties  from  90  to  10^000  kin 

In  this  section  the  method  used  to  derive  a  comprehensive  model  for 
atmospheric  properties  will  be  given  in  detail.  This  model,  covering  the 
altitude  range  from  90  to  10,000  km,  originates  from  the  attempt  to  apply 
the  principles  presented  in  the  preceding  sections  together  with  observational 
data  so  as  to  minimize  the  number  of  assumptions  needed  and  maximize  the 
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accuracy  of  the  calculations.  Major  reliances  has  been  placed  on  using  the 
density  values  from  the  density  model  in  section  h  as  primary  input;  starting 
with  this,  emphasis  has  l)een  placed  on  developing  a  procedure  to  calculate 
the  !-emaining  properties.  The  physical  laws  that  are  assumed  to  apply  are 
the  hydrostatic  equaUon  and  the  equation  of  state.  It  is  also  assumed  that 

diffusive  separation  occurs  above  110  km  and  that  isothemal  conditions  are 
present  above  UOO  km. 

— Derivation  of  500-  to  200-kin  Pressure  Profile  from  Density  Profile. 
The  first  step  involved  in  calculating  the  atmospheric  properties  is  to  derive 
a  pressure  profile  from  the  input  density  profile  (Eq.  2l»)  for  the  200  to 
800  km  altitude  range  and  for  a  given  S'  and  t.  Kallmann  (1959)  has  given 
a  method  for  obUining  a  pressure  profile  from  a  density  profile  in  which 
both  profiles  have  the  same  comparable  accuracy,  although  the  fomer  covers 
a  lesser  altitude  range.  If  (10)  is  integrated  between  two  levels  h.  and 
hg,  where  hg  >  hj^,  the  explicit  relation  for  the  pressure  at  hj^,  P(h,),  is 


P(hj^)  ■  P(h2)  *  j  g(h)^(h)  dh. 


(25) 


where  P(h2)  is  the  pressure  at  hg  and  g(h)  is  given  by  (11). 

From  (12)  and  (13), 

P  •  p  k  T/m  .  (25) 

Using  (26)  and  p  ,  T,  and  M  values  from  Anderson  (1961»),  values  of  P 
were  calculated  near  minimun  (S'  .  70  units,  t  -  5  hr)  and  near  maximum 
(S'  .  200,  t  .  11,  hr)  sunspot  conditions  at  both  000  km  and  SOO  km.  Table  5 
i.  esents  these  values,  together  with  the  ratio  Pgoc/Pjoo.  <*0^  Pgoo  P500 
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•re  the  pressures  at  800  km  and  500  km,  respecUvely.  If  p.  is  taken  equal 
to  sero,  then  the  anall  magnitudes  of  the  ratios  In  this  Uble  Indicate  that 
(25)  can  be  used  to  calculate  a  fairly  accurate  pressure  profile  from  about 
500  to  200  km.  The  largest  error  In  pressure  (6  per  cent)  will  occur  at 
sunspot  maximum  at  500  km. 


TABI£  5e  Pressure  Values  near  Sunspot  Minimian 
(S'  -  70,  t  ■  5  hr)  and  near  Sunspot  Maximun  (S'  ■  200,  t  •  Ik  hr) 


Altitude 

(km) 

Sunspot  Minimun  P 
(dynes/cm2) 

Sunspot  Maximum  P 
(dynes/cm^) 

800 

1.21  X  10“® 

1.35  X  10"^ 

500 

7.03  X  10"^ 

2.3li  X  10"^ 

.02 

^  .  .06 

*^500 

500 

To  evaluate  the  Integral  on  the  right  side  of  (25),  advantage  was  taken 
of  the  fact  that,  for  functions  whose  behavior  Is  approximately  exponential, 
an  accurate  integration  formula  (Kallmann,  195?)  is 


r(h2)  -  F(hj) 

Et'  (hgj  -  In  P  (h^)  ■  *>1)  » 


(27) 
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where  in  this  case  F  ■  g(h)  p  (h)  and  hg  -  h^  is  taken  to  be  a  sufficiently 
anall  altitude  increment.  A  step  size  of  5  km  is  used  for  the  pressure 
calculation  inasmuch  as  computations  made  in  steps  of  1  km  were  not 
significantly  different. 

B.  Initial  Conditions  at  110  kro.  In  this  model  the  altitude  where 
diffusive  equilibrium  begins  is  taken  at  110  km,  thus  the  initial  property 
conditions  must  be  specified  at  this  altitude.  If  the  constituent  concentrations 
and  temperature  are  specified  for  this  altitude,  the  remaining  atmospheric 
properties  can  be  derived  from  them. 

Schaefer  and  Nichols  (1961,  1963),  using  a  massenfilter,  report  values 
of  0.5  at  110  km  for  an  O/Og  ratio;  these  results  are  in  general  agreement 
with  ultraviolet  absorption  results  of  Kiqpperian  et  al  (I960),  according 
to  Schaefer  (1963).  However,  they  are  at  wide  variance  with  previous 
rocket-bome  mass  spectrometer  measurements.  The  major  differences  between 
Schaefer's  and  Nichols'  experiment  and  those  conducted  in  the  past  by  other 
investigators  were:  (a)  the  analyzing  instrument  was  a  massenfilter  instead 
of  a  Bennett  mass  spectrometer;  (b)  the  instrumentation  was  ejected  at 
altitude  from  within  a  volume  pressurized  with  helium  to  minimize  the  effects 
of  rocket-bome  contaminants;  and  (c)  the  ion  source  was  designed  for  direct 
immersion  in  the  ambient  atmosphere  with  the  greatest  possible  open  look 
angle  to  minimize  the  surface  recombination  of  active  constituents.  The 
following  values,  based  upon  the  results  of  Schaefer  and  Nichols  (1963), 
were  adopted  for  the  constituent  number  densities  at  110  km: 
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I 


nCNj)  •  22.0  X  10^^ 
n(0g)  .  5.0  X  10^ 
n(0)  •  2.5  X  10^^ 
n(A)  •  0.16  X  10^^ 

Hence,  the  total  number  denaity  la  n  •  2.97  x  10^^  particlea  cm"^.  The 
density  is  from  (12), 


e 


*^i  ^'i  * 


(28) 


Hhere  is  the  molecular  maas  for  each  of  the  j  constituents.  Using  (28), 
the  density  at  110  km  is  ^  «  1.36  x  10  ^  cm  and  the  mean  molecular 

mass  is  m^Q  -p/n  -  1*.58  x  10‘^^  ga. 

The  mean  molecular  uelght  M  is 


„  28 . 96  m 

M  - - „ 

h.8  X  10 


(29) 


Hence,  subeUtuting  m^^  in  (29),  -  27.73.  The  temperature  at  alUtude 

h  can  be  cos^uted  from  the  molecular-scale  temperature  'L,(h)  at  the  same 

N 

altitude  from 


\  '  (30) 

where  and  are  the  mean  molecular  weights  at  h  and  90  km  altitude, 
respecUvely.  Using  the  value  Tjj(UO)  -  26l  K  from  the  U.  S.  Standard 
Atmosphere  (1962)  in  (30),  Finally,  the  partial  pressure  for 

each  constituent  can  be  calculated  at  110  km  from  (13),  where 

Pi  -  Ui  k  T  .  (31) 
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The  toUl  preaaure  at  110  km  computed  from  (13)  is  p,,.  .  1.02  x  10*^ 

p  XxU 

dynes  cm"  . 

These  initial  conditions  for  properties  at  110  km  will  vary  somewhat 

with  solar  activity  and  diumally.  Inasmuch  as  there  is  a  lack  of  enough 

accurate  measurements  of  these  quantities  at  110  km,  their  variations  cannot 
be  determined  at  present. 

Atmospheric  ProperUes  from  S»0  to  100  km.  Table  6  swinarlzes  the 

above  values  of  the  atmospheric  properUes  derived  for  110  km;  also  presented 

are  the  values  at  90  km  from  the  0.  S.  Standard  Atmosphere  (1962)  used  in 
the  model. 


TABI£  6.  Values  of  Atmo^heric  Properties 
at  90  and  110  km  Used  in  Model 


Altitude 

(km) 


Density  Pressure  Temperature  Mean  Molecular  Number 
(gw  cm-  )  (dynes  cm'^)  (deg  K)  Weight  Density  (cm*^ 


90 

3.17  X  10"^ 

1.61; 

no 

1.36  X  10"^° 

1.02  X  10"^ 

181 

28.96 

6.59  X  10-^^ 

250 

27.73 

2.97  X  10^^ 

The  mean  molecular  wsight  and  temperature  at  any  altitude  h  between  90 

and  110  km  can  be  computed  by  assuming  a  Unear  variation  of  these  quanUUes. 
Hence, 

T  .  T,„  +  ■  ‘‘'90^ 

5>0  - -  ,  (32) 
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and 


•  »*90  * 


(Miio  -  M^q)  (h  -  90) 


20 


(33) 


The  pressure  Pj^  at  altitude  h  is  calculated  from  (10)  and  (26),  thus 

P  m  ff  dh 


dP 


m  ff 

TTT 


Oh) 


Dividing  (3li)  by  P  and  integrating  from  90  to  h  km  gives 


In 


^h  p  mg  dh  ®a  ^  i#  ju  /oe'x 

fZ  •  -s  -h-  --n  f  S'*"*  (35) 

^90  90  ^ 


idiere  g^  is  the  average  acceleration  of  gravity  between  90  and  h  km,  and 
R  is  the  universal  gas  constant,  8.317  x  10^  erg  mole*^  deg“^.  Hence,  P^^ 
can  be  found  by  substituting  Tj^  from  (32)  and  from  (33)  in  (35).  Finally, 
the  density  at  h  can  be  derived  from  (26). 

D.  Atmospheric  Properties  from  110  to  500  km.  Starting  with  the  P^ 
at  no  km,  given  by  (31),  values  can  be  computed  for  other  altitudes  up  to 
500  km  by  (21)  and  (23)  (the  total  pressure  F  is  not  known  accurately  above 
500  km).  However,  since  P  is  known  as  a  function  of  altitude  h  from  200 
to  500  km,  m  (or  M)  can  likewise  be  determined  as  a  function  of  h  for  this 
same  altitude  interval.  This  follows  from  (28)  and  (31),  where 


1  ^  1 
"Al  " 


Pi  m^ 


(36) 


Since  p  ,  m,  and  P  are  known  as  a  function  of  h  from  200  to  500  km,  T 
can  be  calculated  from  (26).  Practicany,  it  is  not  necessary  to  calculate 
T  above  UOO  km,  since  conditions  are  isothermal  above  this  altitude.  Finally, 


n^  can  be  calculated  from  (31)  • 
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However,  p  is  not  known  as  a  function  of  h  between  110  and  200  tan. 
Its  altitude  variation  in  this  altitude  range  is  found  by  first  deriving  M 

as  a  function  of  alUtude  for  this  range.  This  is  done  by  ass««ing  the 
altitude  variation  of  M  to  be  as  follows: 


"  ■  "UO  ■"*(*>-  UO)  +  B  (h  -  110)^  .  (37) 

The  polynomial  (37)  is  a  fairly  accurate  approximation  because  H  decreases 
Slowly  from  110  to  200  tan.  Differentiating  (37)  with  respect  to  h  gives 


(tt 


aj;  -  A  ♦  2B  (h  -  110) 

At  h  -  110  tan,  B  .  0,  and  A  IliJO 

/ 110  — 


-  M 


90 


At  h  .  200  tan.  A  ♦  180  B  \  .  *‘220  *  "200  . 

"an  /  200  -  - -  >  hence 


B~  ^!*220  ~  "200^  -  ("no  -  Hqo) 

3533  - 


(38) 


P  can  be  found  as  a  function  of  h  as  follows:  (a)  starting  with  the 

known  at  110  tan,  solve  (21)  and  (23)  to  find  the  P^  as  a  function  of  P 

from  110  tan  to  200  tanj  (b)  the  m  value  corresponding  to  p  is  found  from  (36) 

(c)  then  the  corresponding  altitude  h  is  detennined  from  (2?)  and  (37). 

To  find  p  as  a  function  of  h  from  110  to  200  tan,  it  is  assumed  that 

P  decreases  exponentially  with  altitude  between  small  altitude  increments, 
thus 


P2  -  Pj^  e‘  “2  (#1  -  *2) 


(39) 
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where  and  Pg  (P^  ^  pressures  at  geopotential  altitudes 

and  respectively.  The  geopotential  altitude  f  is  related  to  the 
geometric  altitude  h  in  meters  as  follows: 

h 

^  J  8  •  (»*0) 

o 


Therefore,  corresponds  to  and  4^  corresponds  to  hg.  From  (39), 


In  (Pg/P^) 

^1  "  ^2  • 


(Ul) 


The  density  pg  at  geopotential  altitude  ^g  is,  from  (10),  (39),  and  (kO), 

P2  ■  *  ■  "‘z  ’’1  *  ■  *2^  •  -  agPj  ,  (1*2) 

Hence,  (la)  and  (h2)  are  used  to  calculate  the  density  from  110  to  200  km. 
The  remaining  property,  the  tenqserature,  can  be  calculated  for  this  same 
altitude  range  from  (26).  The  temperature  value  calculated  at  200  km  can 
be  compared  with  that  previously  derived  to  provide  a  check  on  the  accuracy 
of  the  computations. 

E.  Atmospheric  Properties  from  500  to  8OO  km.  It  was  pointed  out  that 
the  pressure  profile  derived  from  (25)  was  inaccurate  from  8OO  to  500  km. 

The  main  difficulty  was  that  the  pressure  at  the  starting  altitude,  8OO  km, 
was  not  known.  However,  since  the  pressure  at  500  km  is  Known  fairly 
accurately,  (25)  can  be  used  to  calculate  the  pressure  upwards  from  500 
to  800  km.  Also,  since  T  and  p  are  known  from  500  to  8OO  km,  m  can  be 
calculated  from  (26). 
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F.  Neutral  Particle  Concentrations  from  iiOO  to  10,000  Ion.  The  neutral 
atmosphere  above  UOO  km  is  comprised  almost  i*iolly  of  atomic  oxygen,  heliun, 
and  hydrogen;  the  relative  concentrations  of  these  constituents  depend 
strongly  on  altitude  and  temperature.  At  a  certain  altitude  h'  between  500 
and  800  km,  depending  on  temperature,  where  atomic  oxygen  and  helium  comprise 
most  of  the  atmosphere,  the  following  equation  (from  eq.  23  and  31)  holds: 


P' 


**0  -  ^He  *  "He 


k  T  , 


(1*3) 


vrtiere  P*  is  the  total  pressure  previously  calculated  at  h',  and  the  sub¬ 
scripts  0  and  He  are  for  atomic  oxygen  and  helium,  respectively.  Trial 
calculations  indicate  that  h‘  ^600  km  for  S'  <  100  units  (T  1100  K, 

Pig.  U)  and  h'  -700  km  for  S'  >  100  units  (T  >  1100  K).  To  determine  P^ 
at  h'  in  (13),  Uke  (3U)  to  apply  for  only  one  constituent  of  partial  pressure 
P^.  Dividing  (3li)  by  P^  and  integrating  froir  h^  to  hg,  where  h^  >  h^,  and 
using  (11),  gives 


P^Chg) 

TO 


"l  gp  ”o  ^**1  -  _ 

k  4  Ro(hi  +  h2)  ♦  hjhgj 


(Ui) 


Further,  since  the  temperature  and  particle  mass  are  constant  with  altitude. 


^i^^2^  Pi^^2^ 

^i^^l^  ^i  ^“1  ^  ^i^^l^ 


Hence,  P^  is  calculated  from  (Ui),  where  Pi(h2^)  -  Pj^(h')  and  P^(h^)  -  Pq(500  km), 
and  substituted  in  (1*3)  to  find  Pj^^  at  h'.  Using  this  base  value,  Pj^^  and 
rijjg  were  calculated  from  (1*1*)  and  (1*5),  respectively,  from  1*00  to  10,000  km. 
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Alio,  uflng  (U*)  tnd  (US),  Pq  ind  M»re  cilculated  from  500  km  out  to  an 
altitudi,  abovi  about  2000  km,  idiere  atomic  oxygen  becomes  negUgible  in 
concentration. 

According  to  Fig.  1,  the  concentration  of  atomic  hydrogen  in  the 
exosphere  decreases  with  increasing  teii9>erature.  If,  as  a  first  approximation, 
it  is  assumed  that  the  concentration  of  atomic  hydrogen  (atoms  cm"^)  varies 
linearly  with  temperature  T  (deg  Kelvin),  then  at  2000-km  altitude  from  FLg.  1, 

-  5.1  X  10^  -  310  T  .  (1^) 

Using  (Uli),  (1*5),  and  (1*6),  rijj  can  be  calculated  as  a  function  of  T 
from  1*00  to  10,000  km. 

0.  Ion  Concentrations  from  1*00  to  10,000  km.  Although  the  total 
concentration  of  the  neutral  particles  is  at  least  an  order  of  magnitude 
greater  than  that  of  the  ions  out  to  about  2000  km  altitude,  the  ion  concen- 
^***^^°*^®  finally  become  significant  at  the  higher  altitudes.  The  positively- 
charged  part  of  the  ionosphere  above  1*00  km  consists  largely  of  oxygen, 
helium,  and  hydrogen  in  the  singly-ionized  state,  with  o'**  dominant  up  to 
about  800  km,  He”^  dominant  somewhere  between  800  and  11*00  km,  and  dominant 
somewhere  above  11*00  and  1*000  km.  Since  this  ionized  gas  must  be  electrically 
neutral,  each  unit  volume  contains  equal  nianbers  of  negative  charges  (electrons) 
to  match  the  positive  charges.  Therefore, 

"^ion  *  ”e  '  (^7) 

idiere  n^^^  and  n^  are  the  numbers  of  ions  and  electrons  in  the  same  unit 
volune,  respectively. 
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Bowie 8  (1962)  presents  profiles  of  electron  density  over  the  magnetic 

equator  obtained  by  the  incoherent  scattering  technique  from  which  n  can 

0 

be  estimated  as  a  function  of  temperature.  Fig.  5  shows  two  electron 
concentration  profiles  plotted  with  geopotential  altitude,  labeled  1300  K 
and  8uu  K,  derived  from  Bowles'  data.  Inasmuch  as  the  positively  charged 
part  of  the  ionosphere  above  l500  geopotential  km  consists  almost  wholly 
of  H»,  the  temperatures  corresponding  to  the  curves  in  Fig.  5  can  be  calculated 
using  Eq.  ii8,  assuming  thermal  equilibrium  prevails  between  and  the 
electrons.  Thus, 
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g  R 

®o  o 


(hj^  -  hj 
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2  k  In 


V 
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+  R  (h,  + 
0  0  1 
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The  approximate  T  values  of  1300  K  and  800  K  for  the  two  curves  were 
found  from  (1*8)  by  substituting  in  respective  n^^^  (h^)  and  n^^-f  (h^)  values 
from  both  curves  of  Fig,  5,  where  hg  >  h^  >  1^00  geopotential  km,  and,  from 
(1*7),  (h^)  •  n^  (h^)  and  n^^^  (hg)  -  n^  (h^). 

To  estimate  njj+  as  a  function  of  T,  it  is  assumed  that  the  concentration 
of  H-f  varies  linearly  with  temperature  T.  Therefore,  using  concentration 
values  from  the  curves  of  Fig.  5  at  5000  km  altitude  (2800  geopotential 
km  from  eq.  1*0), 

njj+  s  6.0  X  10^  +  8.00  (T  -  800)  at  h  -  5000  km.  (1*9) 


Since  the  positive  ions  below  500  geopotential  km  are  largely  0+,  the  concen¬ 
tration  of  0+  as  a  function  of  T  can  be  derived  in  the  same  manner  as  H+. 
Consequently, 
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-  1.8  X  10^  ♦  920  (T  -  800)  It  h  -  550  km. 


(50) 


To  find  as  a  function  of  T,  concentration  values  were  used  from 
altitude  region  of  the  atmosphere  idiere  the  relative  concentration 
of  He-f  was  greatest.  Accordingly,  njj^-f  values  were  determined  from  the 
two  curves  of  Fig.  5  ^t  1000  km  (865  geopotential  km)  by  graphically  subtracting 
extrapolated  straight-line  portions  of  the  H-f  and  Of  curves.  The  equation 
for  the  linear  variation  of  n^j^f  with  T  is, 

•  1*.0  X  10^  ♦  92.0  (T  -  800)  at  h  -  1000  km.  (Si) 

For  a  given  S»  and  t,  only  one  particular  T  value  is  used  in  (li8) 
through  (5l).  This  is  the  temperature  value  at  1*00  km  that  can  be  taken 
also  to  apply  at  all  higher  altitudes  up  to  about  8500  km.  The  height 
variation  of  an  ionic  constituent  can  be  found  by  (1*8)  in  an  altitude 
region  where  all  other  ionic  constituents  can  be  neglected.  However,  in 
regions  from  1*00  to  10,000  km  where  no  one  ionic  constituent  is  greatly 
predominant,  the  height  variation  of  each  constituent  can  be  found  from 
an  equation  given  by  Hanson  (1962,  equation  8a). 

H.  Atmospheric  Properties  from  8000  to  10,000  km.  Because  the  n^, 
corresponding  to  a  given  T,  are  known  as  a  function  of  h  from  800  to  10,000 
km,  can  be  calculated  from  (31).  Consequently,  P  can  be  found  from  (23). 

Then,  m  can  be  determined  from  (36),  where  n  represents  the  sum  of  the 
neutral-particle  and  ion  concentrations.  Finally,  p  is  computed  from  (26). 
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Th*  equ*Uone  In  SncUon  7  an  pposrMnd  for  an  IM.70»  coapaUr. 
far  IP.  porpoa.  .f  ,Uo.,rau„  .odal,  Uo  .aanplaa  ar.  p™,  pf 

aaaa^  ^  ^ 

»..pa.  -aann.,  m.  o«.r  ...pf,  „„„„„„ 

■ini™.  ».  s.  ,ai„  aaad  r„  ^  p,  ,  ^ 

lar  Paptanbar  fppp,  U»  ....duui,  ,oUo„n,  ».  a„na  M*,.aiU.«d. 
««l.ar  daunauon.  ^  3,  ,3,  _  ^ 

m  and  i^paa.  i»p  (p.  .  f^r  pof,,  3.  .  33  ^  ^  ^ 

iaUoanp  ».  purfdah  Mpn-uut™  aPot.  .oU,  aaanpf.a  a„  f.,  .  .  3.  j,, 

«»  danaiw  for  u,ia  tlna  haa  Paan  „„d  „  oloaal,  app,.aa,al.  ».  dldn^fip 

anrapd  danalty  or  »,  a»  of  tPa  danal.iaa  for  .aa.»  hon,  of  «.  da, 

divided  by  2U. 

The  results  of  the  compuUtlons  made  from  the  model  are  given  in 

two  Ubles  that  present  the  atmospheric  properUes  vs  altitude  from  100  to 

10,000  km.  Tables  7  and  8  exhibit  the  properties  for  S-  -  25o  units  and  S-  . 

I»l»  units,  respectively,  in  addition,  the  pressure  scale  height  H(eq.  llj), 

a  commonly  used  parameter  for  upper-atmosphe«  calculations,  is  shown  in  ^ 
colimn  6. 


yil  ♦  Sunmary 

It  was  pointed  out  in  the  introduction  to  this  report  that  the  approach 
followed  by  most  models  is  to  assuse  altitude  profiles  for  some  of  the 
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TABI£  7>  Upp»r  Ainospbert  NtutriO.  ^pcrtits  Tf«  Altitudi  Ntar  Sunspot 
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properties  In  order  to  derive  the  remainder.  Usually  the  altitude  variation 
of  the  me«i  molecular  mass  was  Introduced  somewhat  arbitrarily^  thus  a 
physically  consistent  vertical  distribution  of  the  composition  could  not 
be  obtained.  To  avoid  this  approach^  the  model  presented  herein  starts 
with  density  profiles  from  an  empirical  density  model  as  primary  input;  the 
remaining  properties  are  calculated  using  knowledge  concerning  the  physics 
of  the  i^per  atmosphere  together  with  boundary  values  derived  from  measurements. 

This  model  cannot  be  described  simply  by  a  collection  of  equations^ 
tables,  and  graphs.  In  order  to  provide  the  proper  basis  for  understanding 
it,  a  description  of  the  upper  atmosphere  Is  given  first  in  Sections  II  and 
III  In  the  form  of  a  concise  discussion  sunmarlzlng  the  present  knowledge 
of  its  composition  and  physics.  Section  IV  explains  the  empirical  density 
model  used  as  a  primary  input  for  the  property  model  in  Section  V.  To 
illustrate  the  model,  Section  VI  gives  the  results  of  computations  made 
with  an  IBM-7090  in  the  form  of  two  tables  presenting  the  average  atmospheric 
properties  for  sunspot  maximum  and  mlnimun  conditions,  respectively.  Also, 
these  tables  can  be  used  to  estimate  the  average  atmospheric  properties  for 
the  months  immediately  following  the  Argus  (September  1958)  and  Starfish 
(July  and  August  1962)  high-altitude  nuclear  detonations.  Tables  derived 
for  many  values  of  S'  and  t  will  be  presented  in  a  later  report. 

Flexibility  has  been  stressed  in  the  constniution  of  the  model  in  order 
to  readily  take  advantage  of  new  information  as  it  becomes  available.  The 
model  is  developed  mainly  for  the  purpose  of  making  predictions  of  atmospheric 
properties  for  comparison  with  measurements  made  from  rockets  and  satellites. 
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Undoubtedly,  the  result  of  such  comparisons  will  suggest  changes  in  the 
model,  hopefully  tending  toward  an  Improved  description  of  the  atmosphere 
wid  more  realistic  model.  Therefore,  the  model  presented  in  this  report 

must  be  considered  the  beginning  of  one  attempt  to  achieve  greater  understanding 
of  the  upper  atmosphere. 
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Appendix  A 

THE  NEUTRAL  DISTRIBUTION  OF  EENSITY  IN  A  MODEL  EXOSPHERE 

Several  authore  have  considered  the  problem  of  extending  the  vertical 
hydrostatic  distribution  of  neutrals  beyomi  the  base  of  the  exosphere 
where  collisions  are  so  infrequent  that  ^Jie  barometric  law  is  no  longer 
strictly  applicable  (Brandt  and  Chambevlain,  I960;  Johnson  and  Fish,  I960; 
Johnson,  1961,  and  Opik  and  Singer,  1959,  1961).  Spherical  sywnetry  is 
assumed  by  each  author  and  collisions  above  the  base  of  the  exosphere  are 
ignored.  Both  Opik  and  Singer  (1961)  and  Johnson  (1961)  classify  the 
exospheric  particles  as  to  orbital  type  and  give  equations  in  the  form  of 
definite  integrals  over  velocity  space  that  describe  their  distribution. 

Their  minerical  results  (obtained  by  numerical  integration)  appear  to  be 
consistent. 

Johnson  divides  the  Maxwellian  velocity  distribution  of  exospheric 
particles  into  four  groups  depending  on  their  trajectories  with  reference 
to  the  base  of  the  exosphere.  These  four  types  are:  (1)  ballistic  elliptic, 
(2)  orbital  elliptic,  (3)  ballistic  hyperbolic,  and  (1*)  orbital  hyperbolic. 
The  corresponding  normalized  groups  are  designated  X2>  Xy  "id  X|^ 
respectively.  The  author  has  found  that  it  is  not  too  difficult  to  express 
these  fractional  groins  in  terms  of  error  functions  and  thus  eliminate  the 
need  for  nimerical  integration.  It  is  the  purpose  of  this  note  to  illustrate 
the  analytic  expressions  for  the  above  distributions  and  to  show  that  quite 
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simple  expressions  can  be  used  to  approximate  exospheric  densities  out  to 
about  5  earth  radii.  It  is  hoped  that  such  a  contribution  will  simplify  the 
use  of  the  results  of  previous  authors. 

It  has  been  argued  that  exospheric  particles  in  group  and  particles 
incident  on  the  earth  in  group  X3  do  not  exist.  One  must  therefore  eliminate 
the  fraction  X|^  +  1/2  x^  of  particles  in  the  hydrostatic  distribution  since 
they  are  not  present  in  the  real  distribution.  However,  there  is  some 
disagreement  as  to  whether  particles  in  the  trapped  orbits  corresponding  to 
group  X2  are  present  or  not.  Brandt  and  Chamberlain  (I960)  argue  that  these 
orbits  are  populated,  while  Opik  and  Singer  (I960)  argue  they  are  not. 

Johnson  (1961)  suggests  that  trapped  orbits  are  probably  populated  out  to 
6  or  7  earth  radii  but  that  at  greater  distances  a  transition  occurs  to  a 
region  free  of  trapped  particles. 

rhe  analytic  e^tpressions  for  the  orbital  types  corresponding  to  the 
definite  integrals  given  by  Johnson  are : 

-  erf  (o(Vg)  -  7T75  (,iv^)  exp  - 

Vl-  ( erf  (« V  Jr7r~J  - 

-v/r/r^)  exp  e*p 

X2  •  vA  -  |erf  (okV^  Vr/r^)  - 

(e(.V^  J^)  exp  (-0k\^/r^)]  exp 

X3  ■  1  -  erf  (oCVg)  ♦  ,4i7TT  (0(V  exp  (-OC^v  2)  . 

"  0 
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^/l  -  -  erf  (ckV^  7r/r^)  -f 

(cs^Vq  ^Av^)  exp  (-<<\^r/r^)^  exp  (-0^^^^), 


^A  -  (r^/r)^  ^1  -  erf  (okV^  ^r/r^)  + 
(dkV^  7r/r^)  exp  (-  Ci^V^^r/r^)\ 


where 


m/2kT, 


Vg^  -  2ga^/r, 

''o^  *  *  '’o^’ 


U 

irf  (Z)  *  VVtt  J*  exp  (-t^)dt, 


r  is  the  distance  from  the  earth's  center,  r^  i^  the  distance  from  the  center 

of  the  earth  to  the  base  of  the  exosphere,  a  is  the  radius  of  the  earth,  g 

is  the  acceleration  of  gravity  at  the  earth's  surface,  m  is  the  mass  of  the 

particle,  T  is  the  temperature,  and  k  is  the  Boltzmann  constant.  The  above 

expressions  can  also  be  expressed  in  terms  of  the  scale  height  H,  since 

«.^V  ^  -  r/H. 
e 

The  following  combinations  are  also  useful: 

Xi  +  X2  •  erf  («<v  )  exp  (-ot^  ^), 


XB  +  Xj,  '  1  -  erf  (ol  V^)  *  .-CT"  WVg)  exp 
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X2  *  Xjj  -  «A  -  exp 


(7) 


and 

(8) 

As  one  might  expect,  X2  X|^  ®re  independent  of  since  these 

two  groups  comprise  the  total  fraction  of  particles  with  velocities  less 
than  escape  velocity  and  greater  than  escape  velocity  respectively,  and  this 
property  is  clearly  independent  of  a  reference  sphere.  The  purpose  of  a 
reference  sphere  is  to  permit  distinction  between  orbital  and  ballistic 
particles.  It  is  interesting  to  note  that  the  use  of  a  reference  sphere 
introduces  a  discontinuity  at  the  base  of  the  exosphere  since  from  (3) 

1/2  X^(r^)  >  0.  However,  this  discontinuity  is  very  small  and  indicates 
that  ignoring  collisions  above  this  somewhat  arbitrary  sphere  is  a  very  good 
approximation.  The  above  expressions  are  very  slowly  varying  functions  of 

r  . 

0 

With  trapped  orbits  fully  populated,  the  fraction  X^^  *'■  X2 
of  the  hydrostatic  distribution  represents  the  actual  distribution.  The 
fraction  x^  1/2  x^  of  the  hydrostatic  distribution  represents  the  ballistic 
distribution  (trapped  orbits  absent).  It  can  be  readily  seen  from  Eqs,  (1), 

(2)  and  (3)  that  the  analytic  expressions  for  the  fractions  X^^  X2  1/2  X3 
and  +  1/2  X3  somewhat  unwieldy.  Herring  and  Kyle  (I96I)  and  Aamodt 
and  Case  (1962)  have  previously  given  the  lengthy  analytic  expressions  for 
the  special  case  of  a  ballistic  density  distribution  in  terns  of  error 

a 


Xi  +  X3  “  1  -  exp  . 
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functions.  However,  it  turns  out  that  is  much  smaller  than  Xi  below 
about  5  earth  radii.  Therefore,  below  this  altitude  it  is  pe:iiissible  to 
neglect  1/2  snd  to  use  the  much  simpler  expression  (5)  and  (8)  to  correct 
the  hydrostatic  distribution. 

The  hydrostatic  density  distribution  is  given  by 

f  (r)  •  f  (Tq)  exp  ^  (1  -  '•/r^)J  . 

We  have  then  approximately 


51  [i- 


TTFTtT^  exp  ^  (r) 


-  ■\4I7nJXVg)  exp  ^  (r) 


where  ^  is  the  density  assuming  only  ballistic  orbits  are  present  and 
p^(r)  is  the  density  with  the  trapped  orbits  fully  populated 
The  overestimate  in  using  (10)  and  the  underestimate  in  using  (11)  amounts 
to  only  about  10^  and  h^.  respectively  at  a  distance  of  5  earth  radii  from 
the  center  of  the  earth  and  with  an  isothermal  temperature  of  12^0  K.  The 
error  increases  slowly  with  increasing  temperature  and  distance  from  the 
earth.  It  is  not  too  difficult  to  justify  the  use  of  (10)  and  (11)  since 
(a)  they  are  very  good  approximations  below  5  earth  radii,  (b)  it  is  unlikely 
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that  the  trapped  orbits  are  either  completely  populated  or  completely 
unpopulated  (the  above  approximations  give  results  that  lie  between 

wd  and  (e)  the  expressions  (1),  (2),  (3)  and  (1*)  are  only  approximations 

based  on  a  simplified  model. 
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Fig.  1 

Fig.  2 

Fig.  3 

Fig.  1* 
Fig.  5 


.  Distribution  of  atomic  hydrogen  vs.  altitude  above  the  earth  at  the 
extremes  of  the  sunspot  cycle  (Johnson,  1961a).  The  curve  labeled 
1^00  K  is  for  the  average  temperature  corresponding  to  the  maximum 
of  the  sunspot  cycle,  while  the  curve  labeled  1000  K  corresponds 
to  the  sunspot  minimum. 

.  Fraction  of  Maxwellian  particles  in  the  atmosphere  vs.  altitude 
for  three  temperatures  and  two  cases. 

.  Penetration  of  solar  radiation  into  the  atmosphere  according  to 
Friedman  (1959).  The  graph  shows  the  altitude  at  which  a  fraction 
e"^  remains  of  the  radiation  incident  on  the  atmosphere. 

.  Isothermal  temperature  (above  UOO  km)  vs.  S»  for  every  hour. 

.  Electron  concentrations  vs.  geopotential  altitude  for  temperatures 
of  750  K  and  1300  K,  respectively  (Bowles,  1962). 
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